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ABSTRACT 



^ ' There are now a half dozen young pulsars detected in high energy photons by 



the Compton GRO, showing a variety of emission efficiencies and pulse profiles. We 
present here a calculation of the pattern of high energy emission on the sky in a model 
which posits 7-ray production by charge depleted gaps in the outer magnetosphere. 
This model accounts for the radio to 7-ray pulse offsets of the known pulsars, as well 



compute the polarization angle variation and find that the previously misunderstood 



> 

l/^ I as the shape of the high energy pulse profiles. We also show that ~ 1/3 of emitting 

^ ■ young radio pulsars will not be detected due to beaming effects, while ~ 2.5 x the 
number of radio-selected 7-ray pulsars will be viewed only high energies. Finally we 

o 

^\ . optical polarization sweep of the Crab pulsar arises naturally in this picture. These 

(— I , results strongly support an outer-magnetosphere location for the 7— ray emission. 

o 



Subject headings: pulsars — gamma rays — polarization 
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^ • 1. Introduction 



Since the birth of 7-ray astronomy the two most prominent galactic point sources have 
been identified with young rotation-powered pulsars, the Crab and Vela pulsars. The site of 
7-ray production has however been intensely debated for almost two decades. The two principal 
models are the 'polar cap' picture exemplified by the work of Daugherty and Harding (1982) and 
the 'outer gap' model championed by Cheng, Ho and Ruderman (1988, CHR). In the former 
acceleration occurs near the neutron star surface and the observed emission results from a 7 — S 
pair cascade. In the second model acceleration is posited in charge depleted regions near the 
speed of light cylinder and 7 — 7 pair production is an important process. We have examined 
these two models (Chiang and Romani 1992) and find that the observed pulse profiles arise most 
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naturally in a modified version of the outer gap picture. Moreover, spectral calculations (Chiang 
and Romani 1994, CR94) have shown that in this model emission processes vary throughout the 
magnetosphere and can produce spectral variations through the pulse like those seen. In this 
paper we quantitatively compare the outer gap predictions with the observed profiles, supporting 
these conclusions. 

The impetus for these computations is the dramatic increase in information on the 7-ray 
pulsars provided by the Compton Gamma Ray Observatory. In addition to improved light curves 
and phase resolved spectra for Crab and Vela {e.g. Nolan et al. 1993), CGRO has detected at least 
four other pulsars in high-energy emission, including Geminga, which has not been seen in radio 
emission {e.g. Thompson, et al. 1992, Ulmer, et al. 1993, Fierro, et al. 1993, Mayer-Hasselwander 
et al. 1994). Other bright plane sources await confirmation as pulsars, and the expected number 
of each of these classes of objects is presently unknown. It is apparent from these data that 
older rotation-powered pulsars are increasingly efficient producers of 7 emission, at least for 
characteristic ages ^ 10^'^ — lO^y. In deciphering the puzzle of 7 production efficiencies upper 
limits on other pulsars are also important. Their interpretation is, however, unclear. Finally 
detailed information on the location of the emission region is an essential tool for understanding 
the polarization and spectrum seen from the observed pulsars. The geometrical computations of 
pulsar beaming in this paper address each of these issues. 



2. Emission Zone Geometry and Profile Calculation 

Well above the neutron star surface the dipole component of the magnetic field will dominate. 
We have computed the location of the outer magnetosphere charge depletion regions (outer gaps) 
by calculating the dipole field structure from the full retarded potentials. The outer gaps (CHR) 
are regions above the null charge surface (where the co-rotation charge density changes sign) 
extending to the speed of light cylinder, where co-rotation must break down. Accordingly plasma 
lost from this gap cannot be replenished from the stellar surface and a charge depleted region with 
significant E • B is developed. This is the location of the particle acceleration which gives rise 
to the observed 7-ray emission. Following CHR, we assume that the outer gap will be bounded 
below by the surface of last closed field lines. We determine this surface by finding the dipole field 
lines which are tangent to speed of light cylinder. Traced back to the stellar surface, these give a 
polar cap moderately extended in longitude for significant magnetic inclinations a. The upper gap 
surface will be along field lines for which substantial 7 — 7 pair production shorts out the outer 
gap. It is this upper surface, where the population of accelerated pairs builds up to appreciable 
densities, that in fact constitutes the emission zone of the gap. 

This emission surface lies significantly within the open field line region above the polar cap. 
We find {c.f. Chiang and Romani 1994) the soft photons supporting the 7 — 7 pair production that 



limits gap growth are provided principally from regions in the inner magnetosphere (the lowest 
portions of the gap, where synchrotron losses are significant) or from the stellar sm'face. Pair 
production will be strongest close to this soft photon source, and so we assume that the gap width 
for a given position around the polar cap is determined by this pair production at the null charge 
surface. Two effects determine the pair production length and hence the gap width at this point. 
First as one moves to the sides of the polar cap, the intersection of the null charge surface with the 
last closed line moves higher in the magnetosphere. There the soft photon flux at the threshold for 
pair creation on the curvature 7 from primaries accelerated to radiation reaction limited energies 
in the gap will decrease - this tends to increase the pair production length. Countering this to 
some extent, the angle between the central soft flux and the locally produced high energy 7's will 
increase higher in the magnetosphere - this lowers the pair creation threshold and for a steep 
soft spectrum from thermal or synchrotron radiation, decreases the pair production length. The 
complete solution of the gap growth and radiation production is a complex non-linear problem 
(CR94). Here we adopt a simple scaling law for the gap width w at the null charge surface which 
approximates these behaviors w oc re, with Tc the local field line radius of curvature and with 
the width at the gap's lowest point scaling as a~^. This lowest width can also depend on the P 
and B, thereby controlling the pulsar luminosity - this has not been treated in detail. Thus gap 
regions starting at high altitude will be wider. However, the gap width does not grow as fast as 
the B rs^ r^'^ spreading of the field lines; the flux threaded towards the sides of the gap decreases. 
Along any given field line, charges will flow above the null charge surface, so we take the set of 
field lines determined above to define the gap upper surface (Figure |l|). 

We note that we have not included the effects of particle inertia on the field lines or the field 
perturbations induced by magnetospheric currents. These should be small in our near vacuum 
region as long as we are not very close to the light cylinder. In particular these produce only small 
changes in the geometry of the polar cap. Moreover, we wish to emphasize that the majority of 
our results are quite robust to perturbations of the field geometry and are thus independent of the 
details of how we have chosen our gap surface. However, fine structure in observed pulse profiles 
can help detail the precise location of the high energy emission zone. The brightness along the gap 
builds up as the particle density increases, but then falls off as the light cylinder is approached. 
As an estimate of this variation we take F oc 2s — s'^ , where s is the arc length along the field line 
measured in units of the light cylinder radius, and impose an exponential cutoff at s = 1 with 
a = 0.5. 

With the emission zone defined, we note the radiating particles will have significant 7 factors 
along the field lines. High energy photons will thus be radiated in a narrow beam tangent to 
the emission zone. As E^ decreases below an MeV, an increasing fraction of the emission will 
arise from lower energy pairs and the inverse-Compton and synchrotron photons producing the 
observed pulse will have a larger spread about the field line (CR94). However, significant fiux will 
still be along the local field line and so the high energy beam shape is defined by the tangent to 
this last closed surface. Aberration effects and light travel delays will be very significant high in 





Fig. 1. — Two views of the surface of last closed field lines, cut at 1.2rLc for an a = 70° rotator 
(see text). The gap lies above the null charge surface, inside the open field line region. We show 
the illuminated section of the the gap upper surface. The observed pulsar 7-rays are beamed along 
this surface. The neutron star's radius is increased to O.Ivlc for clarity. 



the magnetosphere, so we account for these and bin the observed radiation into pulse profiles. We 
display this as a skymap in pulsar phase and observer co-latitude (C = a + /?)• In Figure ^ we 
show the radio caps (from both hemispheres) along with the 7-ray beam from the Southern radio 
cap. Angles and the indicated line of sight are those appropriate to the Vela pulsar. 



2.1. Gamma Ray Beam Characteristics and Radio Data 

Several features of the 7-ray beams are worth mentioning. First, since radiation occurs 
above the null charge surface, where fi • B = 0, 7-ray emission observed at a given (^ comes from 
the polar cap in the opposite hemisphere. Particles of the opposite sign from those producing 
the observed emission will be accelerated inward in the gap. Their radiation would be visible 
in the same hemisphere as the radio cap, however CR94 show that the radiation efficiency and 
spectrum can be quite different for the inward-going particle; we expect that they will produce 
little observable flux. 

Second, as noted by Chiang and Romani (1992), double pulses with significant bridge emission 
are visible for most C ^-s long as a is not too small. These pulses are caused by crossing caustic 
surfaces in the pulse phase sky map, where flux from a large portion of the outer gap arrives in 
phase. These caustics are produced only with the full effect of aberration and light travel time 
and can be quite sharp when the emission surface is thin. Relatively narrow 7-ray pulses will be 
visible for lines of sight close to the upper edge of the gamma ray beam, and will have two poorly 
defined peaks with strong bridge emission, such as PSR1706-44 and possibly PSR1055-52. Also 
very wide (width between 7-ray peaks A > 0.5) double pulses are possible in this model, especially 
for small a and C ~ 90°. As a approaches 0°, regions where the outer gap approaches the light 
cylinder after crossing over the rotation pole become increasingly important in the pulse profile. 
Faint emission from this region is visible on the trailing side of the pulse in Figure ^. Near a = 0° , 
the 7 emission is concentrated to the equator, and the region where strongly pulsed radiation is 
visible will be small. Those pulsars seen at small a will tend to have wide pulses. 

Radio data are particularly important in constraining possible viewing geometries. In 
particular, polarization data can give very good determinations of /? from the rotating vector 
model (Lyne and Manchester 1988, LM). Pulse widths can also be used to estimate a in some 
cases (Rankin 1990, 1993), but in general a is only poorly constrained, unless emission is visible 
in an interpulse or through a large range of pulsar phase. However, the radio data (especially with 
accurate polarization measurements) also give good measurements of the phase of the magnetic 
axis. The relative phase of the radio and 7 pulse are now available for the observed pulsars (c./. 
Ulmer 1994). In Figure y we compare the 7-ray pulse width (between peaks) as predicted from our 
model as a function of the phase offset from the radio with the values measured for the observed 
7-ray pulsars (Table 1). The agreement is very good, although in the case of the Crab, the radio 
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Fig. 2. — Skymap of pulsed beams and pulse profiles for Vela parameters (a = 65°). Above: Pulsar 
emission as a function of pulsar phase - observed lines-of-sight cut horizontally across the image at 
constant (; the Earth line-of-sight is indicated (dashed line). Emission from the observed outer gap 
(halftone) and both radio caps (contours) are shown. For clarity, the outer gap visible for (" > 90° 
is not shown. Below: Pulse profiles for Vela. The model 7-ray pulse for the line of sight above is 
indicated, along with the leading radio pulse and schematic plots of the hard (magnetospheric) and 
soft (thermal) X-ray pulses. 
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Fig. 3. — T-ray pulse width A vs. lag of first 7 pulse from the radio pulse 5. Observed values 
for radio pulsars are indicated. Model predictions for a covering the observed pulsar's range are 
given by the full lines. The phase of the magnetic axis for the Crab is somewhat uncertain - the 
horizontal bar gives 1/2 the range over which the core component is expected to appear. 



emission lies slightly closer to the 7 pulse than expected. However, we note that for this shortest 
period pulsar, the envelope in which the radio emission, identified with the precursor, can appear 
is relatively large. Half this range is shown by the horizontal bar. For the other pulsars more 
complex profiles and polarization data allow the true centroid of the radio pulse to be better 
determined. Overall, the good reproduction of the observed correlation is a strong success of the 
model. 



2.2. Individual Pulsars 

We have computed models for a range of a and sampled the pulse profiles of these for the 
full range of (,. Figures ^a and |jb show contours of constant b and A respectively in the a — C, 
plane. The gap width used in these sums is appropriate to a ~ 0.1s pulsar, near the middle of the 
observed range. Above the line of A = 0, no 7-ray emission is seen. Radio emission is visible for 
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Fig. 4. — Isocontours of lag of from the radio phase b (left) and 7-ray pulse width A (right) in 
the a — C plane (see text). Contours are computed for a gap width appropriate to a median age 
(P r^ 0.1s) pulsar; selected contour values are labeled near a = 90°. Observed pulsars are positioned 
at estimated values of (a, C); for Geminga a large range along A = 0.5 is allowed. Pulsars known 
to have (^ > 90° (f) have been displayed reflected to the Northern hemisphere. Above the line A = 
no 7 rays are visible; far from /3 = pulsars will not be detected in the radio. 



small (3; lines showing the range of visibility for the Crab are drawn, although at low frequencies 
and/or low intensities, radio flux may be visible over a wider range. Interpulsars, with a near 
90° and radio emission visible from both poles are in the lower right corner of the diagram. For 
a > 90° the geometry will be similar with /? of oppposite sign. 

Using /3 values (and a, where available) from the radio data we have plotted positions of 
young pulsars on Figure HK and Hb. In some cases, we have chosen a value of a consistent with 
the radio and 7-ray data. For Geminga a substantial range along A = 0.5 is consistent with the 
data. Results for each of the pulsars are in good agreement with the models. Moreover several 
peculiarities of the observed profiles are explained with these figures. For example, it is apparent 
that for much of the parameter space giving A ^ 0.5, as for Geminga, /3 will be too large for the 
radio pulse to be visible. Also for nearly aligned pulsars such as PSR1929+10 (Phillips 1990), no 7 
pulse is expected, so the strong EGRET upper limits can not be construed to imply a luminosity 
cut-off at the characteristic age of PSR1929-I-10. Also, in several cases pulsars are located just 
outside the range of 7 visibility {eg. PSR1509-58 and possibly PSR0354-I-55); in these cases, the 
GeV 7-ray beam tightly collimated to the gap surface may be missed, while the lower energy 
emission with its wider beam pattern may be detected. In these cases OSSE, or even X-ray 
energies, may provide the best hope of pulse detection. 

Most importantly, this model gives significant predictions for the behavior of the radio 
polarization. Careful measurements should be able to check these expected viewing geometries 
in many cases. Conversely, polarization studies in connection with our high energy model 
computations can help decide when careful pulse searches of young pulsars are likely to be 
successful. Ultimately, knowledge of the viewing geometry will help in locating the emission region 
for various sectors of the pulse profile and, through phase-resolved spectra, will allow a detailed 
study of the emission region. 



2.3. X-ray and Optical Pulses 

In this paper we do not compute X-ray (0.1-10 KeV) profiles in detail. However, we note that 
there are several important sources of X-ray emission in this model. First we expect the region of 
the gap near the null charge line to be a significant source of X-ray flux for young pulsars. This 
flux will radiate widely and will be visible as a broad 'hard' pulse. We predict that such flux 
will appear at phases ~ 0.35 — 0.5 after the radio pulse. Ogelman (1994) has noted such a hard 
'magnetospheric' pulse component in young pulsars and has shown that it connects spectrally 
onto the high-energy power-law emission. A second X-ray pulse will result from softer flux with 
a roughly blackbody spectrum emitted from hotter regions near the magnetic axis on the young 
star's surface, due both to anisotropic magnetic opacities and to heating by backflowing particles 
from acceleration zones (Arons 1983). These soft pulses will lie near phase 0, but should again 



10 



be quite broad, due to wide emission zones and general relativistic bending of the photons near 
the neutron star surface. The radio emission at ~ 10 — 100 stehar radii should not suffer this 
effect. Several young pulsars, have in fact been observed to have hard and soft X-ray pulses 
shifted in phase, in broad agreement with this picture {eg. PSR1055-52: Ogelman and Finley 
1993; PSR0656+14: Cordova et al. 1991, Finley et al. 1992; Geminga: Halpern and Ruderman 
1993; and Vela: Ogelman 1994). When a is close to 90° emission may be visible from both caps 
and from both gaps, complicating the profiles. A third, primarily unpulsed soft X-ray component 
can be present due to the initial cooling of the hot neutron star; following Ogelman (1994) we 
expect this to be significant only for young (lO^y < Tc < lO^y) neutron stars where the surface T 
is still large, but the gap luminosity is relatively small. 

At present optical pulse profiles are available only for the Crab, Vela and PSR0540-69 pulsars. 
For the Crab, correspondence of the optical and 7-ray profiles indicates that optical emission is 
produced over much of the outer gap. However, for the Vela profile the optical pulse is substantially 
narrower than that of the high energy emission. The phase and width of the pulse are consistent 
with an origin near the null charge surface (c./. Figure |2|). PSR0540-69 has a broad pulse with 
two peaks (Middleditch and Pennypacker 1985, Gouiffes et al. 1992) consistent with a line of sight 
cutting fairly high on an outer gap. Optical detections have however also been established for 
Geminga (Halpern and Tytler 1988) and suggested for PSR0656-I-14 and PSR1509-58 (Caraveo et 
al. 1994a, b). For the former two pulsars the optical emission may be from the Rayleigh- Jeans tail 
of the thermal flux from the heated polar cap, although non-thermal gap emission is preferred. 
PSR1509-58, however, if correctly identified at rrir ~ 21 must be dominated by non-thermal gap 
emission. This counterpart is bright enough to allow determination of an accurate pulse profile; in 
our model we would expect a single broad peak of emission aligned with the hard X-ray pulse {i.e. 
at phase ~ 0.3 with respect to the radio peak). 



3. Crab Polarization Results 

Interpretation of the Crab pulsed emission has long proved problematic. It is clear from 
the alignment across a wide spectrum that the radio-optical-X-7 emission from the main pulse 
and interpulse arises from similar zones, which we identify with the outer magnetosphere gap 
(although the intensity weighting along the gap surface will differ significantly for the different 
energy bands). We follow other authors in identifying the precursor with the normal radio 
pulse. Confusion arises however because the optical pulse exhibits a strong double sweep in the 
polarization position angle - this has long been interpreted as surface emission in a two-pole RVM 
fit with a ~ 90° {e.g. Narayan and Vivekanand 1982). Since the Crab is the only source for which 
polarization information is presently available from the high energy region, this association has 
been very influential in the interpretation of other 7-ray pulsars. 

Adopting the same premise as the rotating vector model, (Radhakrishnan and Cooke 1969) 
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namely that the polarization vector of the emission is set by the local magnetic field, we can 
similarly compute the position angle swing through the pulse in our outer gap picture. We 
compute the vector radius of curvature of the field line at the emission point, project this on the 
plane of the sky, follow aberration variations, and average the position angles from all regions 
contributing to a given pulse phase to obtain the model polarization sweep. Remarkably, because 
of the combination of emission from a range of altitudes in the magnetosphere, our model produces 
an initial sweep from the leading edge of the outer gap (the first pulse), a reversal and a final 
sweep in the second pulse. This behavior is obtained precisely for the viewing angles that produce 
the observed ~ 140° profile width. The left panels of figure || shows the optical polarization 
position angle variation (from INT data. Smith, et al. 1988) along with a high energy pulse profile 
from BATSE. Figure ^ also shows model results with the viewing angle chosen only to match 
the observed pulse width and separation. While it is gratifying that the pulse profile matches 
well to the data, it is remarkable that the corresponding polarization sweep shows a striking 
similarity to the observed data, as well. We consider this a major success of our model. With 10m 
class telescopes and modern detectors it should be possible to obtain polarization data for Vela, 
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Fig. 5. — Crab pulsar profiles and polarization sweep. Left: light curve data from BATSE 
(Fishman, GRO Newsl.,1, 6), INT optical polarization data (Smith et al. 1988). Right: Model 
results for Crab outer gap. 
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PSR1509-58 and PSR0540-69, as well. These would be particularly important in motivating an 
emission location in the lower portions of the outer gap - since double sweeps generally occur only 
for a ^ 70°, C ^ 70°, we expect most of these pulsars to show a single sweep of position angle. 



4. Population Estimates 

In addition to allowing acceptable models for individual pulsars, it is important for a 7-ray 
emission picture to reproduce the population as a whole. To test our model on this front, we 
have assumed that magnetic inclinations are random at birth and that the observer orientation 
is also random and have integrated the beaming fractions and detection probabilities for the 
radio and 7-ray emission over a population of young (r < 10^'^y) pulsars. For the radio emission 
we have adopted conventional beaming fraction relations {W = 6.5°P^^''^ LM; W = 5.8°P^^'^ 
Rankin 1993). The results for the fraction of all young pulsars detectable in the two bands are 
shown in Figure ^ for the LM beaming law; values for Rankin's beaming law differ by only a few 
percent. We have also considered how the results might be affected by alignment torques secularly 
decreasing a, but with typical alignment timescales ~ lO^y (LM) the differences in detection 
fractions are very small. Uncertainty in the P — P dependence of the gap width makes these 
fractions slightly uncertain; improved estimates of gap closure can remedy this. 

Our population estimates show good agreement with present limited data. First, as is well 
known, the mean radio beaming fraction of these young pulsars, 0.27, leads to a pulsar birthrate 
reasonably consistent with the galactic supernova rate. Second, roughly 1/3 of the radio selected 
pulsars should not be visible in 7-rays, to a consistent luminosity threshold. Of the young radio 
pulsars monitored, several do in fact have upper limits on their efficiency of 7-ray production 
lower than that of the known emitters (Thompson 1994). Thus with ~ 5 radio pulsars detected in 
7-rays it is not surprising that significant upper limits exist for PSR0656-I-14 and PSR1929-I-10; 
non-detection for ^ 2 pulsars should be due to 7 beaming. There is also a very strong upper limit 
for PSR1951-I-32, but this is probably a mildly recycled pulsar and has a dipole field roughly an 
order of magnitude lower than the young pulsars. Because of decreased synchrotron emissivity 
and difficulty in gap closure and inverse Compton scattering photon production, we would not 
expect this to be an efficient 7-ray pulsar. True millisecond pulsars are similarly poor targets for 
magnetospheric 7-ray emission. 

Finally, because of the large 7-ray beams, the existence of objects such as Geminga not 
detected in the radio is not surprising. Indeed, scaling from the the 5 radio pulsar detections, we 
expect ~ 12 non-radio pulsars to be visible in 7-rays to flux levels comparable to the faintest radio 
selected object, in agreement with the observationally motivated estimates of Helfand (1994). We 
thus expect most of the ~ 20 unidentified galactic point sources to be young pulsars. In most 
cases identification of the pulse period will be very difficult: X-ray emission, such as that used to 
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find the period of Geminga (Halpern and Holt 1992) offers the best hope, but most sources will be 
3-10 times more distant, so the count rates will be low. Support for the pulsar hypothesis may be 
obtained by finding the expected hard Inverse Compton Scattering spectra above lOOMeV, or in 
some cases by association with old supernova remnants. 

It is important to note that theories placing the principal 7— ray production in acceleration 
gaps at the polar cap (Daugherty and Harding 1982, Harding, Ozernoy and Usov 1994; Sturner 
and Dermer 1994) have great difficulty in meeting the population test. These models cannot 
easily produce the range of pulses with bridge emission observed. In particular, to produce wide 
pulses these models require very small q ^ 10° with the observer line of sight C nearly aligned 
with the rotation axis ft to intercept this small pole (Harding and Daugherty 1993, Dermer and 
Sturner 1994). Even if pulsars are preferentially born with aligned magnetic axes (a hypothesis 
not supported by the radio data), the present significant number of detections coupled with the 
small beaming fraction in these models already implies an enormous birthrate of young pulsars, 
typically ^ 30x the galactic supernova rate. We consider this to be unacceptable. 

Thus while some 7-rays may be produced in polar gaps, this cannot account for the emission 
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Fig. 6. — Venn diagram for pulsar detectability (beaming fractions). Random magnetic 
inclinations, high energy emission from both poles and radio beaming factors from LM are assumed. 
Values are from an integration over pulsars with Log(-B) = 12.5 and Tc < 10^'^y. 
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presently observed from radio pulsars. The energetic arguments used in these models to estimate 
7 fluxes will carry over to the outer gaps in many cases. The processes converting the primary 
current to the observed 7 emission will however be quite difl'erent. Accurate models are not yet 
available for the outer magnetosphere (CR94), so that luminosity laws for 7-ray pulsars must at 
present remain phenomonological. 



5. Conclusions 

We have computed the expected pulse shapes and beaming fractions for 7-ray emission from 
young radio pulsars in a model of charge acceleration in the outer magnetosphere. The results bear 
a strong resemblance to 7-ray pulse profiles detected by CGRO. In addition our model reproduces 
the phase relationships between the radio emission and the 7-ray pulse. Observed X-ray and 
optical pulsations find a natural interpretation in our geometric picture, as well. Important 
connections between radio pulse and polarization properties resulting from the viewing angle and 
expected 7-ray profiles have also been established - these give useful predictions for individual 
pulsars and allow strong tests of the model. We further find that the polarization properties of 
the Crab pulsar, with emission arising in the upper magnetosphere are easily explained in our 
model; this supercedes previously confusing interpretations of the polarization data. Finally, the 
numbers of objects detected in the radio and 7-ray channels are as expected for a population of 
young pulsars and we infer that most of the unidentified galactic plane sources will be pulsars. 

These results are based principally on the geometry of the emission region and are thus 
relatively free from uncertainties in details of the emission process. We are in the process of 
generating quantitative models for the radiation produced in the upper magnetosphere, but 
complicated non-linear models will be needed to give accurate results. Previous estimates of outer 
gap spectra and fiuxes (CIIR88, Ho 1990) give rough estimates of the total energy available, but do 
not give accurate results for the observable flux or spectrum. Nonetheless, the dominant physical 
processes in this outer gap have already been identifled by CIIR88. The substantial amendments 
needed for a realistic model, however, mean that spectral results are not yet available. 

We feel that the success of our geometrical sums flrmly establishes the location of the 7-ray 
production in the outer magnetosphere and thus resolves the long standing debate with polar cap 
models in favor of an outer gap picture. Unless similar results can be duplicated by polar models, 
we feel that these models are not viable and efforts to compute the spectrum and luminosity from 
the polar cap site cannot explain the bulk of the observed pulsar emission. Much work remains 
to be done to provide a full picture of the origin of pulsar gamma rays, but assignment of the 
radiation site gives good hope for further progress. In particular CR94 show that the radiation 
processes vary strongly with altitude, and that these variations can be mapped directly to pulse 
phase. Thus phase resolved spectra provide a keen diagnostic for refining global emission models. 
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Ultimately, detailed comparison with observed profiles should allow us to probe the inertia and 
current perturbations that are not followed in the present calculation. Understanding of these 
will help greatly in unraveling the mechanics of the pulsar acceleration process and in producing 
a self-consistent model of the pulsar magnetosphere. Thus 7-ray measurements can become an 
important tool in deciphering the puzzle of the pulsar phenomenon. 

Support for this work was provided by NASA grants NAGW-2963 and NAG5-2037. We are 
grateful to Jim Chiang and Joe Fierro for important discussions on the EGRET observations that 
improved this paper. We also thank F. Graham Smith and Jamie Biggs for sharing the Crab 
optical data. 



Table 1: Young Pulsar Parameters 
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PSR 



P(ms) log B{G) logT(y) logE'(erg/s) 



L^/E" 



A" a{°f /3(°) X-rays<= 



7-ray Pulsars 



B0531+21 


33 


12.6 


3.10 


38.65 


0.004 


0.05 


0.40 


80 


-18 


h 


B1509-58 


150 


13.2 


3.19 


37.25 


0.012-'' 


0.31 


~0 


60 


-15^^ 


h 


B0833-45 


89 


12.5 


4.05 


36.84 


0.012 


0.11 


0.43 


65 


+14 


s,h 


B 1706 -44 


103 


12.5 


4.24 


36.53 


0.024 


0.21 


0.32 


42 


+25 


s?,h 


B0630+17 


237 


12.2 


5.53 


34.54 


~0.1-0.3 




0.49 


25^=^ 


+50^^ 


s,h 


B1055-52 


197 


12.0 


5.73 


34.48 


0.2 


0.22 


0.32 


70 


-8 


s,h 



Candidates 



B0540-69 


50 


B1951+32 


40 


B0656+14 


385 


B0355+54 


156 


B1929+10 


227 


B0950+08 


253 



2.7 


3.22 


1.7 


5.03 


2.7 


5.05 


1.9 


5.75 


1.7 


6.49 


1.4 


7.24 



38.17 
36.57 
34.58 
34.66 
33.59 
32.75 



.33 



8 +8 

51 +4 

31 +20 

170 +5 



s.-^ 



"assumed beaming fraction 0.5 

''radio pulse - first 7 ray peak separation (phase units) 

^^7 pulse width (phase units) 

''estimated from 7-ray data, when available. Values uncertain for Geminga. 

'^s=soft flux from surface/heated polar cap, h=hard flux from base of outer gap 

■'^0.1-1 MeV flux only - not detected by EGRET 

^estimated from X-ray, optical, radio data. 
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